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ABSTRACT 


The Generalized Rotor Performance (SRP) program isa 
computer program designed for calculating forward flight 
performance of a helicopter rotor system at a specific 
flight condition. It can be used tc evaluate either an 
articulated or a hingeless single rotor system in forward 
flight or in a wind-tunnel test. Тһе program was originally 
designed by the Sikorsky Aircraft Company and purchased by 
the United States Navy. 

The goals of this thesis were (1) to reinvestigate the 
theory and logic used in the program, (2) to add selected 
desirable features to the program, (3) to produce a much 
needed Users' Manual, and (4) to run an analysis comparing 
the program's calculated results against manufacturer's 
data. These goals were accomplished and the results of the 
analysis indicated that the program produces highly accurate 
results within the normal cruise range of a modern 


helicopter. 
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I. INTRODUCTION 

The Generalized Rotor Performance (GRP) program is a 
computer program designed for calculating forward flight 
performance of a helicopter rotor system at a specific 
flight condition. It can be used to evaluate either an 
articulated or a hingeless single rotor system in forward 
flight or in a wind-tunnel test. The GRP series of programs 
were ocigi nally designed by the Sikorsky 20 ٣ 
Corporation with the cooperation of the United Aircraft 
Research Laboratories. Work began ОП this series ОБ 
programs in estando has continued to date. The United 
States Navy purchased this program in 1964 and has used it 
in the Naval Air Systems Command as a helicopter performance 
computer program. 

The goals of this thesis were (1) to reinvestigate the 
theory and logic used in the program, (2) to ensure that the 
Navy's version of this program performed the calculations 
according to the design theory, (3) to produce a much needed 
Users' Manual, (4) to add certain desirable features to and 
correct errors found in the program, (5) to run an analysis 
comparing the program's calculated results against 
Manufacturer's data and (6) to document areas that may need 
further attention. The work was done with the help and 
cooperation of both faculty at the Naval Postgraduate School 
and personnel in the Naval Air System Command. 

The output available from the program includes rotor 
SHAE: НОСОМ, rotor profile horsepower, main rotor 
Bere ዓን drag, chrust, H force, rolling and pitching 
moments and other forces and moments calculated in the 
Shaft, control and relative wind reference axis systems. 
Also available are the azimuthal histories of the blade's 


flapping angles, rates and accelerations, plus azinuthal and 


oN 


DUL 


шиг СОг ес ОЕ angles of attack, CL, ср, Mach number, 
sweep angles and air load distribution. In addition, the 
program will calculate a Fourier coefficient series for 
flapping angle, radial station air lcads and azimuthal 7 
Force distribution. 

The program can be divided figuratively into three main 
routines. The first routine determines a steady state 
flapping solution. The second routine determines the forces 
and moments generated by this flapping soluticn. The third 
routine compares the calculated results of lift, drag and 
other options to those desired by the user. pf Zeche 
calculated values are not within a predeternined tolerance, 
the program will, by use of a first order Taylor series, 
generate new values to reenter into the blade flapping 
routine. 

The program uses a combined blade element/momentum strip 
theory. The blade's radius is divided into a maximum of 15 
segments. The blade is advanced around the azimuth ina 
prescribed number of degrees. At each particular azimuthal 
position, velocities are computed in the perpendicular (UP), 
tangential (UT) and radial (UR) directions. From these 
velocities and a known pitch angle distribution, the lccal 
angles of attack are calculated. From these angles of 
attack, the forces at each radial blade segment аге 
determined. Once the forces are known, the local moments 
can be found. From the summation of these moments at a 
particular azimuthal роза топ, the blade flapping 
acceleration is calculated. This acceleration, combined 
With the calculated flapping angle and rate, is used to 
advance the blade to the next azimuthal position. A steady 
State flapping condition is assumed tc exist when the blade 
flapoing angle and rate at the PSI equal zero and 360 degree 
azimuthal position are within a prescribed tolerance of each 


other. 





Pars nethodmor Calculation accounts for retreating blade 
stall, the reverse flow region and compressibility effects. 
The program can accommodate a full range of geometric and 
design variables including flapping hinge offset, elastic 
flapping hinge restraint, first and second harmonic cyclic 
inputs and spanwise variations in blade twist, local nass 
densities, chord and tip sweep. The program uses no small 
angle assumptions and has no restrictions on tif speed, 
forward velocity or advance ratio. The rotor system can be 
oriented in any direction in space andcan be given any 
rotor shaft or Meera E TOll, pitch ог yaw angular 
velocities. The program will perform calculaticns in 
straight and level flight or a uniform induced velocity nay 
be added to simulate climbs and descents. The GRP will 
accept up to five different airfoil data decks plus one 
blade spar characteristics data deck. LIFE 6 
information is entered into the program in the form of CL 
mee versus Angle or Attack. Tables for up to 15 different 
Mach numbers. The user has a choice of six methcds of 
soluticn depending upon the desired restrictions. The user 
has available nine printout options, two of which аге 
program debugging options. There are also error printout 
messages that will assist the user having difficulty with 
the progran. 

In the method described above, the flapping motion 
determined is about a flapping hinge offset, and with the 
elimination of all assumptions that the flapping angles are 
small, the rotor control axis (axis of no feathering) and 
ше сір path plane axis (axis of no first harmonic flapping) 
are no longer considered convenient for reference in the 
analysis of rotor blade потоп. Therefore, the axis 
selected for use in the analysis is the rotor shaft axis 
system. All forces, moments and angıiss are referred to the 
Shaft axis system with the exception of some of the final 
output forces which are referred to the relative wind and 


control axis systems. 


ld 








With the use of the computer, the GRP program eliminated 
many of the simpliflying assunpticns of the earlier 
classical theory originated by Wheatley and Bailey in Ref. 1 
and 2. The assumptions of the classical theory did not 
impose serious limitations in low speed flight, but as 
helicopter sveeds increase, the inaccuracies inherent in the 
theory seriously limit the usefulness of this method. The 
assumptions of the Wheatley-Bailey theory which are not 
present in the GRP program include the following. 

1. The flapping and inflow angles are assuned 

small. 

፡ Рре lift and drag coefficients are approximated 
by a linear and a quadratic variation, 
respectively, with angle of attack. 

Meee The effects of Mach núnber on CLE and CD are not 
considered. 

4. Sectional characteristics in the reverse flow 
region are the sane as those in the 
conventional flow. 

5. The sectional characteristics of blade twist, 
tip sweep, flapping hinge offset and root cut 
Out are ignored. 

While the GRP program represents a refined approach to 
the classical theory, there are still some basic assumptions 
which make the GRP subject to error. These are: 

1. Steady state two-dimensional airfoil data are 

used. 

2. Quasi-static blade analysis is used. 

3. The rotational speed about the shaft axis is 
constant. There 15 по lead cr lag action in 
the progran. 

4. The rotor blade is assumed rigid in bending and 


torsion. 
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5. The rotor inflow is assumed uniform unless 

varied by the user. 

6. Spanwise flow 15 псог рогатей into the 

calculation of blade angle of attack; however, 
it is assumed that the flow at one segment does 
not affect the flow at any other segment. 

It is felt that the errors induced by the above 
assumptions are relatively small in the normal cruise speed 
region of a modern helicopter. This is the region from just 
below the minimum power airspeed to the maximum allowable 
cruise speed. The program can not calculate rotor hover 
power. This is because one of the factors in the 
denominator in the main routine which estimates new flapring 
routine reentry parameters is the advance ratio. Since the 
advance ratio is zero in a hover, the computer would attempt 
to divide here by the number zero. Highly accurate results 
are not available in the airspeed region from hover to just 
below the minimum power airspeed while using the normal 
uniform inflow assumption. This is а region of highly 
Hen-Uniform ficw. Variable inflow may be inputed by the 
user. However, at this time, no information is available on 
how successful this technigue is in accurately eStimating 
the required rotor system horsepower. While the program's 
required technique for inducing harmonic variable flow is 
cumbersome, Bramwell in Chapter Fcur of Ref. 6 describes 
several methods which could be incorporated into the 
progran. 

ከከ ጠው በጋን) Static blade analysis aSSumption is valid 
except in the very high speed region where a large 
percentage of the retreating blade is in a stalled 
condition. This is an area of aerodynamic hysteresis that 
is influenced by unsteady aerodynamics. While encountering 
lift hysteresis, the amount of lift change above the steady 
State CL is about the same as below. However, the lift 
distribution on the rotor disc will change. Seal. „us 


delayed and occurs at a Slightly later azimuthal station 


13. 
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рдеп ртейістей by quasi-static analysis. smie difficult, 
complicated and computer-time-consuming procedures can be 
taken to reduce these assumptions there is no guarantee that 
the accuracy of the solution will improve. It is felt that 
the present program produces highly accurate results in the 


flight range of interest in a modern helicopter. 





A. GENERAL COMMENTS 


The GRP uses a combined blade element/momentum strip 
theory in its calcultions. The blade radius is divided into 
a maximum of 15 segments or strips. Figure 1 illustrates a 
typical blade element. In crder to cbtain a steady state 
flapping solution the program must calculate the time 
history of the flapping motion. This necessitates a 
complete knowledge of the flapping angles, rates and 
accelerations. The time history cf this motion must be 
found by solving the highly non-linear eguation for the 
Summation of moments about the flapping hinge. Figure 2 
illustrates the forces which produce moments on the rotor 
system. The blade has aercdynamic forces acting upwards, 
blade weight acting downwards, centrifugal ferce acting 
Ошечавйз, an Optional elastic hinge restraint acting to 
return the blade to a zero flapping angle and an inertia 
force resisting any change in blade flapping. The moment 
Summation equation about the flapping hinge may be written 


as 


The moment due to inertia force can be expressed as re, 
where I, is the blade moment of inertia about the flapping 
ПЕЕ cubstituting chis ito equation 1 and solving for 0, 


the following governing equation is obtained. 


(2) o = Hs&&o + Un t Mee + Tene 
I, 


кыс аа OE attack and local Mach number at sach 


blade segment were known, the CL and CD could be obtained 


I 








= ሂ/ዬ 
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Figure 1 Blade Element Diagram 
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Figure 2 Rotor Blade Mcment Diagram 
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Ул 


from the inputed airfoil inicrmation. The local aerodynanic 
forces and moments can be calculated and the flapfing 


acceleration determined. Equation 2 can be rewritten as 


K seqme nis 


(3) Фа = ጄ ( амь: + Мы; + амр. + аМеца: ) 
421 Ig 


The flapping angle and rate at the (N + 1) azimuthal 


position could be obtained fron the following expressions. 


(4) ort ES En Ы Фе at 
° 225 as 
(5) Ban = Bn + nat ት Ва > 


However, since the flapping 1S periodic in nature and 
has a direct relationship to the azimuthal angle, PSI, the 
values for flapping are solved with respects tc PSI, vice 
time. The values of £1, v and time are related by the 
equation AW =naT. Therefore 


Е x 
(6) б = 88 . nd = TL E 
At ач 44 
ee > de 
(7) В = à? = пав = -2- 5 
ó ++ ay” 
The governing equation for the flapping motion now 
becomes 
$ XK 
(8) б = Magen t+ ቁ. + Mer t Nene 


2. 
en 


The flapping solution is based on the assumpticn that 
the angle of attack is known. However, it is not and the 
program must proceed through an iterativ2 process in corder 
len DD CAC inflow Tatio, collective pitch and cyclic 


input angles required to generate the desired forces. 


i. 


В. ANGLE OF ATTACK CALCULATIONS 


A very important and basic part of this program is the 
procedure by which the local angles of attack ате 
calculated. While the program will calculate angle of 
attack with any angular velocity applied either to the rotor 
system or the helicopter, the development here will describe 
level flight only. The classical approach ignores radial 
flow, UR, and the angle of attack would te obtained as shown 
in Figure 1. However, as the blade rotates abcut the shaft, 
it will encounter a large variation in radial flow. The 
program attempts to compensate for this radial flow in the 
following manner. Instead of the inflow angle PHI equalling 
the arctangent of UP/UT, it is set equal to UP/UTUR wheres 
UTUR is the resultant velocity in the tangential and radial 
direction. СЛЕЗЕ Таа геа іп Figure 3. The pitch 
angle is also reduced by the cosine cf the sweep angle. The 
angle of attack is now calculated in the sweep plane. This 
three-dimensional angle of attack is lower than the 
classical two-dimensional angle. 

The program enters the CL, CD tables with this sweer 
plane angle of attack and the sweep plane resultant Mach 
number. The program computes the forces using the velocites 
in the sweep plane, UP and UTUR, and the blade chord 
geometry in the normal plane. Once the forces are computed 
in the sweep plane they are resolved into their respective 
directional forces. 

This three-dimensional angle of attack, due to sweep, 
Will delay stall on the rotor by reducing the angle of 
attack. This describes what actually occurs on the blade. 
However, it is felt by previous personnel using the progran 
that there is a point where the sweep becomes so large that 
it tends to wash out the lift being produced. The program 
has been modified to reduce CL by one half for sweep angles 
between 60 to 372.5 degrees and reduce CL to zero for sweep 
angles between 72.5 and 90 degrees. Theses high sweep angles 


occur normally only оп the inboard tlade segments of the 
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Figure 3 Spanwise Flow Diagram 
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retreating blade in and near the reverse flow region. 
ОР, UT, UR and pitch angle for level flight are shown in 
Figures 4 and 5. They are calculated in the following 


Manner. The radial velocity, UR, is calculated as 
(9) UR = (Усоза) созЧ созв 


The tangential velocity, UT, has two components. The 
first is the local rotational velocity, Ar, the second is a 
sinusodial component of the forward flight velocity. The 


general expressicn for UT is 
(10) UT -.fLr * (Vcosdg)sin' p 


This expression contains a small angle assumpticn in the 
Borse for blade flapping angle. The program acccunts for 
the fact that the flapping angle does reduce the true radius 


slightly by using the following formula 
(11) ОТ = (E/R + (r - E/R)cosB)N + (Vcosas)sinW 


The perpendicular velocity, UP, consists of three terns, 
the inflow ratio, a flapping velocity and a small component 


purorwdrd velocity. The inflow ratio is defined as 


(12) A = Msinds - v 
R 


The seccnd component is a verical flapping velocity 
Ween 15 2 LUNnet20n Of flapping rate and radius. This is 


computed as 


(13) ПР (2) = (г - Е/В) В 
The third component is due to the fact that there isa 


small component of the radial flow which acts in the UP 


direction due to blade flap angle. This is equal to 


а. 
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9 
}— elk —}+—___( R- E/R) cos Q — 


Cb) 


ри ( V cos Xs) Cos iy 


( N cosas 
Ё (3 соѕ ds) SIN Y 


UT = (E/R + (г ~ Е/Е)созб)а + (Усова) зіп Ч 


Figure 4 UT Diagram 
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UR = (Усов 4%) созфсоз() 
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Figure 5 UP and UR Diagram 
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(14) ПР (3) = (Tcosds)cosY 511 
ine total formula for UP is 
(15) UP = ኣ.ሴ ፪ር656 мн — =/ጾ)ሮ + (Vcosds) cos¥sin® 


The pitch E (0) is expressed in equation 16. 0.75 
is the pitch angle at the 75 percent radius position, 8' is 
the twist depending on the relationship of the location of 
ээ олаје segment to the r = .75R location. The next four 
terms are the first and second harmonics of cyclic pitch and 
(tan $ )@ ls the coupled effect of tne flapping angle on 
the pitch angle. 


(16) 9 = (845 + @ - AlScos¥ - 21Ssin¥Y - A2Scos2W 


- B2ssin2y - кап (6) 0) COS 7 


The local Mach number is calculated as U/a where а 15 


the speed of sound and U is given in equation 17. 
| 
(17) U = (UP> + UT? + ur)” 


The angle or attack can now be calculated. Figure 1 
illustrates that 4 <- ©-ф . Initially, the program requires 
estimated values for the inflow ratio, Collective pitch at 
шэг аата! zero position, harmonic cyclic inputs, blade 
Most and Initial flapping angle and rate. The user can 
either input these values or accept the program's automatic 
default values Of -.02, 5, 1.2, ШІ ت80‎ 0 
respectively. Using these assummed values, the initial 
angle of attack can be determined as fcllows. 

UP = AL Rcos® SCE = 5/፳)6 + (Я со5 45) cosYsinf 
reduces to 
UP = ALR 
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ጠበ Е a/R) cOspjt. + (Vcosds)sinY 
reduces to 
UT 


ПЕ ИШЕ S E/R) N LL 


UR 

reduces to 
UR = Vcos es 

In the program Vcosd, is replaced by the tern „aR where u 


(Vcos ds) costcos$ 


is the advance ratio in the shaft axis system. Since the 
local pitch distribution and inflow ratio are estimated, the 
local angles of attack can be determined. The next section 
describes the method used for calculating flapping angle and 


rates at the N + 1 azimuthal position. 


-> መን ፡፣ፐጥ ከኮ ርንቨንንም SOLUTION FOR THE FLAPPING EQUATION 


In the preceding section, equation 8 was developed in 


Seder tO Calculate the time history of the flapping» moticns 


ችት 
(8) ( = Maeno * Mu * Mee + Mene 
TS 
The relation involves a complicated second order 


differential equation for establishing the flapping angle as 
a function of PSI. The numberical solution is accomplished 
by use of a finite difference equaticn and a step-by-step 
procedure. An important characteristic o£ the solution is 
chat 1% eis periodic in nature. The ВАСЕ 1 0000000 — 
represents a steady state flapping sclution has the property 
<). Using this fact, a Pourier harmonic series 


CN EFE یپ‎ ተራስ to describe the blade flapping mction. 


(18) (2 = AO - Alcos¥ - Bisin¥ - A2cos2y - E2sin2Y - 
со ህወ Во... 


Ву asSuming that the first harmonic flapping is much larger 
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than the other higher harmonics, the series can be reduced 
го 

(19) 6 = AO - Alcos¥ - 815119 

This equation can be differentiated with respect to PSI 
to obtain ۹ 

(20) 0 AlsinY - Bicosy 

(21) СМ = Alcosw - Bisind . e" 

Assuming that the values of 6, @ and 6 are known at 


some azimuthal position, the following equations must hold 


ዘ 


(22) Bun = እ0 ፦ А1соз9 - Bl1sinw 
% 

(23) 9% = Alsin» - В1со5% 

(24) ұм = Alcos¥Y + Blsin% 


These equations can be solved for the N + 1 azimuthal 
position by substituting Way, = ዥ ሌሎ into the above 


formulas. The flapping angle equation becomes 
БЕП 88 Ево а1соз(Щ FAY) - Elsin(%, +AQ) 


By using the following two identities equation (25) can 


be rewritten as eguation (28). 


(26) cos (W, +A) COSY, COSAY ~- sink sin ay 


(27) ኋከ (፡%። там) 5104 соБАЧ + cosh Sinay 
(28) Ena = A0 - A1(cosMcosay -sint sinay) 
- Bi(sinY,\cosAy - cos sinAy) 


The terms can be rearranged into equation 29. The same 


+= 
procedure can be used to develop equation 30 for В. 


АО - cosA¥(Alcos¥, + 81510.) 
= + sinaw(Alsin¥n - Bilcosy) 
(30) A = cos AW(Alsin'Pa - B1lcosWa) 
+ Sinay(Alcos + 31sinw) 


(29) mar 


Substitution into equations 22, 23 and 24 reduces these two 
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expressions to the flapping equations used in the progran. 
The user can either enter the value for AW or accept the 
programs automatic default value of 15 degrees. 
که‎ xx 
(39 Ana = бо со5 49 + Ва Sin Ач 
2 T 
(32) (5 ሓኒ = да ~ Cx sw Av * G- Cos AN) En 


While this integration scheme is nct one of the standard 
methods used, it is very useful in obtaining periodic 
solutions for differential equations similar to the one used 
here. Notice that for small values of AW the trigncnetric 
expression can be reduced to an ordinary Taylor series. By 
assuming sinAW equals AY and cosAW equals one equation 29 
reduces to 


Зе A ہت‎ 
(33) 520 = ба ~ AY on 


By assuming sin Aw equals ÀY and ccsAÀWM equals the first 
ሬ 
twc terms of the cosine series 1 - ሌጁ equation 32 can be 


reduced to 


4 ۱ , 
G9 Qu, - 6ዕ ጐ ሌኛ ክሬ * m S 


ШІ FORCES, MOMENTS AND RADIAL INTEGRATION 


The forces acting on a rotor blade may be found by the 
Summation of the elementry forces along the span at any 
aZimuthal position. The forces considered by the progran 
are the resulting aerodynamic forces only, and are initially 
Summed in the shaft reference axis system. The program 
computes forces in three axis systems. They are the (1) 
shaft, (2) control and (3) relative wind axis systems. 

INE BESTE cadjally integrates differently for lift and 
drag calculations. The drag is calculated from the hince 


онов, Б/В, to the rotor tip. Пади iS calculated iron 
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the hinge offset to the next to last rotor blade segment. 
The last segment is considered a blade "Tip Loss Factor" 
segment. It is assumed that the tip trailing edge vortices 
cause no lift to be produced in this segment. The normal 
procedure is to define this segment as the last three 
percent of the rotor blade radius. 

The first of the maximum 15 blade segments is 
considered the spar or cut out segment. This is defined as 
the area between the hinge cffset and the point where the 
airfoil actually begins. If no spar data are entered, it is 
asssumed that this first segment preduces no lift and drag 
1s obtained by using the CD verse Alpha Tables for the first 
blade airfoil data deck entered. If soar data are entered, 
the first segment and all other segments designated spar 
segments will have the lift and drag characteristics of the 


entered spar data. 


Е.  AZIMUTHAL INTEGRATION METHOD 


Once the integration of the rotor forces and moments 
along the blade is complete, an integration around the 
azinuth must be performed in order to obtain the average 
forces and moments. Since the solutlen of the flappang 
equation is obtained by a step- by- step method, the 
integrands of the integrals over Y аге known only at a 
certain number of equally spaced points around the azimuth. 
Beorecqurlibrium flapping the integral is a periodic function 
of Y , and for this case integration can be shown tc be 
equivalent to an averaging process. This result makes the 
azimuthal integration simple. The following method is used 
where N is the number of azimuthal positions used in the 


Ssalewlations and b is the number of blades. 


ጌኛ 8 


мож 
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КОРА ۹۹0 ITERATION 


Once the program has calculated a steady state flapping 
solution and has determined the resulting forces and 
moments, a question remains to be answered. Is this the 
desired solution and, if not, what must be done to obtain 
Enos sclution? In order to solve the flapping equations 
certain known or assumed values were used, including inflow 
ratio (A), collective pitch (9.75) and the cyclic pitch 
IS, BIS, А25 апа 825). One method, first described in 
Ref. 3, is to iterate on the required lift and drag thrcugh 
modification of A and 9.75, where 8.75 is the pitch angle 
at the 75 percent radius station at the PSI equal zero 
azimuthal position. The modified values are then reentered 
into the flapping routine and the calculation is repeated 
until it converges to within a specified tolerance on the 
required lift and drag. Drag is used here in the sense of 
negative rotor propulsive force. The program procedure is 
outlined below. 

mie required rotor lift and propulsive force are 
expressed in terms of the magnitude and direction of the 
resultant force in the longitudinal plane. These are shown 


in Figure 6 where it can be seen that 


(36) а= о a' 


+ Ц 
(35 вр= ( (О. + (а) »" 


(38) L = R cos (a ") D = Ry sin (a") 
The required lift and propulsive force and their 


resultant are Shown in Figure 7. In a Similar way, 


(2) R xm + Da ) 


(40) ас arctan (Da /La) 


28. 


The differences between the required and the computed Re 


and a" values are defined as 


| 


(41) A 54 a 53 Ry 


(42) ሌ 8"= an E 


to compensate for the 


Ti order to correct A and 8.75 


difference between (Re. ale) and (Rp o") , the required 
кес аге expanded in a Taylor series with A and 9.75 as 


variables. Тһе first order equations are: 


(43) Rog = Ro + DR se 
М 

(44) Qe = ој + YA AA + за! лос 

ХА Bar 


tor the iteration on እ and 6.75 


Solving the equations 


yields 
32" ARg - Me да" 

(45) КА - 3875 9 Өлз 
de за? 2 342 Эс 
3^ 


ФА зе-< جو یج‎ 


385 ла" — da ልዊ 


(56) S Bas => EN А 
ود‎ ја" _ 384 Эс 
ом 39645 3604s 95 


" and a" are known, the corrected 


Now that 92.” Ro and ato 
values of A and 9.75 can be approximated by: 


(47) Пас እኤእ Ваза Pus = 565 


the values of the 


In order to solve equaticns 45 and 46, 
found. The 


partial derivatives in these equations must be 


шо. 
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Figure 6 Resultant GRP Force Calculation Diagram 





Figure Y Required GRP Force Diagram 
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procedure used in based on the Wheatley-Bailey method and 
the formulas can be found on pages 186 and 207 of Ref. 4. 
Reference 5 outlines the derivation and a complete 
derivation was performed and verified in conjunction with 


this thesis. 


G. FLOW CHART 


The flow chart in Figure 8 is a block diagram showing 


the relationship between the various parts of the GRP 


program. 
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111. GRP_USERS' MANUAL 
The GRP currently has the capability to enter over 200 
individual case variables and option selectors. Presently, 
170 are available to the users. This manual explains the 
ከ ከ ከ ከ ከ and Output format of this GRP progran. Examples of 
both input and output are given in order to make the use of 
this program easier. This manual is divided into the 
following areas of discussion. 
A. Main Iteration Options 
В. GRP Data Deck Order 
с. Airfoil Data Deck Requirements 
D. Spar Data Deck Requirerents 
E. Sample Data Input Format 
Е. Case Input Listings 
G. Case Input Default Values 
ДС сасе Input Data 
I. Case Input Fornat 
J. Case Optional Output Indicators 
Ка Во 360 Execution Control Cards 
L. Sample Program Output 
It is recommended that the user examine the GRP Case 
Input listing carefully, since a few options require certain 
variables to be inputed which are not necessarily lccated in 
the same general area of the input listings. An attempt has 
teen made to list all input variables concerned with the 
different options in the discussion of each option and input 


variable. 


ИЕЛІК hep RATION OPTIONS 


The GRP offers the user six different options for 


determining a solution. They are as fcllows. 
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1. Normal Routine 


— — — —— ہے‎ A: ASA A s Qs жы 


In the normal solution to the problem the computer 
Weer Vary the inflow ratio, LAMBDA, and the pitch angle, 
БИ шп огдег to produce the required lift and drag. The 
ас inputs, AS and 815, аге considered constants. 
Uniform inflow is assumed unless the user induces a 
non-uniform inflow by the use of variakles 117 and 118, LAML 
and UVL. The program will calculate the required shaft axis 
angle and position of the control axis from the fixed value 
SE BIS. In all of the options, the user has a choice of 
using either a program calculated first estimate of flapping 
angle, velocity and acceleration or an initial value of zero 
for all of the above. Either method usually requires the 
same approximate amount of computer time. The variable 
PCNV, item 97, determines the method tc te used. 


2. Desired Flapping Angles 


This option allows the user te specify the desired 
longitundinal and lateral flapping angles with respect to 
the rotor shaft axis. The program at each intermedient 
mee iteration will vary A1S, B1S, inflow. ratio and pitch 
angle. Variables 100 and 112, 113 and 114% control the use 


s his option. 
3. Short Iteration Scheme 


This option follows the normal routine with one 
Major exception. The program will make only one pass 
низового спе flapping routine on each major iteration. The 
program may or may not arrive at a steady state solution for 
flapping in the first few iterations. It is assumed that 
the first few iterations in the normal routine are cnly 
rough estimates of the way the variables should be changed 
and that an exact flapping sclution is not actually required 


if the user is only interested in transient flapping 
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behavior. In order to use this routine, the user must input 
a negative number for the variable  XITLIM, item 73. The 
absolute value of XITLIM will still determine the maximum 
allowable number of times the program will enter the major 


iteration routine. 
Ч. Trimmed Moments 


This option will follow the Ncrmal methcd described 
in paragraph one. The Normal method „will only converge on 
lift and drag and will not consider the moments produced. 
If the variable TRIM, item 160, were assigned a non-zero 
value, the program would attempt to trim out the pitching 
and rolling moments about the rotor shaft. It is suggested 
that this be done in a two case Tun. The first case should 
be a normal run, with the desired printout. Then, ror the 
Second case, input the variable TRIM. This will do two 
things. First, it will provide a ccnverged solution without 
consideration of moments. Secondly, it will allow the full 
number of iterations to be used +o reduce the moments. Tne 
program does this by a short routine varying А15 and BIS. 
During this process, the whcle flapping and force iterations 
must be repeated, but it will at least start by using a 
converged solution. The variable PCNV, item 97, should be 
non-zero to allow the flapping solution from the previous 
case to be used as a first estimate in this case. Setting 
Mo variable SKIPIN, tem 91, equal to zero will alla 
force and moment summation to be outputed fcr each major 
iteration. This will allow the user to see exactly how the 


program is proceeding. 
5. TOP Option 


In this option, the program iterates ирсп the 
required lift but ignores any values inputed for required 
drag. This option can be used to simulate a wind-tunnel 


test. The user must input the shaft angle, item 111, and 


CO 
سا‎ ° 


he pitch angle, item 87. These two inputs will be held 
constant. The program will iterate on the inflow ratio, 
LAMBDA, in order to obtain a solution. The variable ТОР, 


ከ ከ ከወን controls the use of this opticn. 


6. ALOPT Option 


— 


This option, like TOP, is a wind tunnel option. It 
also iterates on required lift only. However, here the 
shaft angle, item 111, and the inflow ratio, item 88, are 
required inputs and are held as constants. The program will 
iterate on the pitch angle, item 87, in order to determine a 
Semution. The variable ALOPT, item 110, controls the use of 


BEES option. 


B. GRP DATA DECK ORDER 


Data is entered into this program in the following 
order. 
Peev LLI Coefficient Table 
22 11220131 Drag Coefficient Table 
(Repeat steps one and two as necessary.) 

3. Case Input Data 

ንይ ۶7071669۶ ط۲‎ 67160 Ratio ** 

=e Spar Lift Coefficient Table жж 

6. Spar Drag Coefficient Table — 


** Optional Data 


ئ٠ ت۸30 ع۶۶۷7‎ DECK REQUIREMENTS 


The GRP program requires that all CL and CD information 
ше “entered into The program in tabular form. Tables for up 
to 15 different Mach numbers and five different airfoils can 
be entered. The program currently does not use cr require 
values for the Coefficient of Moment, CM. Since certain 
segments of the retreating blade are in the reverse flow 


région, angle of attack tables are required to include 
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valued from -180 degrees to +180 degrees. If they are not 
included, an error message will be printed when the progran 
can not locate a value for CL and CD at these large positive 
and negative angles of attack. If complete angle cf attack 
information is not available for a particular airfoil, the 
user can use the values provided in Section E. It is 
realized that available data on airfoil behavior at large 
angles of attack are very limited, but so is the region on 
the rotor disc where the blade operates at these high 
angles. Since this occurs only immediately arcund and 
Within the reverse flow region where dynamic pressure is 
low, little precision is lost in performance calculation by 
using one common representation for most airfcil behavior. 

AS a minimum, two values for CL and CD at two different 
angles of attack and Mach numbers must be supplied. As a 
maximum, 15 different values of Mach number may be entered. 
Each Mach number may contain up to 48 different values of CL 
and 48 different values of CD and their associated angles of 
attack. 

The first Mach number must be equal to zero. This table 
can be an exact duplicate cf the lowest Mach number table 
the user has available. The highest Mach number table 
should be high enough in order to prevent the program from 
stooping because of a local Mach number higher than that in 
the table. A quick check can be obtained by adding together 
the rotor tip velocity and the forward flight speed. This 
combined velocity, divided by the local speed of sound, must 
be less than the maximum Mach number entered into the 
program. The program linearly interpolates between Mach 
numbers and angles of attack in order to determine the value 
Of СІ апа СИ. The subroutine BLINY does the intergolaticn. 

Several options are available to help reduce tne number 
of data points that must be entered. If the airfoil is 
Symmetrical, the user only needs to enter values for 


positive angles of attack. The proaran will assign the 
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ammroplate sign to the value of CL and CD according to the 
Sign of the angle of attack. This is accomplished by case 
input variable 107, SYM. This option can also be used for a 
cambered airfoil where values of CL and CD at negative 
angles of attack are unknown. 

Values for large angles of attack need not be entered 
for each Mach number by making use of the program's input 
variables 156 and 157, or their automatic default values. 
Values for large positive and negative angles of attack need 
only be entered for the two lowest Mach number tables. The 
lowest Mach number table must be at a Mach number equal to 
zero. If anangle of attack is greater than variable 156 
(HIALFA) or lower than variable 157 (LOALFA), or the 
programs default values of plus and minus 30 degrees, the 
Mach number is set equal to zero. This ensures that cnly 
the first two Mach numbers have to carry the whcle range of 
angles of attack from -180 to +180 degrees for a cambered 
airfoil or from zero to +180 degrees for a symmetrical 
твої і. 

The format of the table input will now be described. 
The first data card contains the variable 1۰ This 
controls whether or not the user receives an echo printout 
of the Blade Airfoil Data being entered. If tne value of 
WBLADE is equal to zero, the user will not receive an echo 
printout of the Blade Data. If WBLADE is a ncn-zero number, 
the user will receive the echo printcut. The read format 
шог (2 ГпА0Ё 15 F10.0.  WBLADE is the only item оп the first 
data card. 

The next card also contains only one piece of data, 
ከከ NBLADE is the number of different airfoil data sets 
to be used and appears on this card in an 12 format. This 
program will accept up to five different blade airfoil data 
sets. It will also accept one blade spar data set. IE the 


rotor blade being analyzed were composed of three different 


іт 


types of airfoils, NBLADE wculd equal three. However, i 


38% 


саас consisted of three sections, of which the first 
and third section were the same, NBLADE would equal two. I> 
is explained later how the blade segments are assigned their 
respective airfoil type. This arrangement provides the user 
with the ability to vary the make-up of the blade while only 
having to enter into the program once a particular set of 
airfoil data. 

The above two variables, WBLADE and NBLADE, are only 
entered once for each complete computer run which uses the 
same set of airfoil data. Тһе following information will be 
entered twice for each type of airfoil used. It will be 
егей first дог CL and secondly for CD for each type 


TOLL. The overall format is summarized below. 


Card 1 WBLADE 
Card 2 NBLADE 
Card 3- CL's for airfoil number one 


Србешешна ст еј | number ons 
Cis ROrralle cre number two 
CD's for airfcil number two 
Repeat as necessary. 

Card number three contains the variable NZ, whicn is the 
number of Mach numbers for which CL's will Бе entered for 
A erst airfoil. This number must Бе right justified in 
12 format. The maximum number of Mach numbers for each CL 
ир tor one airfoil is 15. This is the only пипрег 
entered on this card. 

Card number four begins the actual Mach number, CL 
versus angle of attack tables. The format in this paragraph 
must be repeated for each Mach number. This first card 1s 
даи ие слее а t¡Welds (12, 10F7.0). The fest field is а 
ОН, ШЕ gustrfled integer in 12 fofhat. It is 
equal to twice the number of data pairs for tnis Mach number 
plus tyo. This tells the computer how many numbers are 
required to be entered for this particular Mach number. А 


data pair consist of one angle of attack and its asscciated 
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лог Ср. The remaining ten fields are ach seven columns 
Hong an floating point or F/.0 format. These fields begin 
Мишо епа У. 10, 17, 24, 31, 38, 45, 52, 59, and 66. On 
BS first card, the field that starts in column three 
contains the number of data pairs at this Mach number. The 
field that begins in columns 10 is the actual Mach number. 
The remaining eight fields on this card are for the first 
four data pairs starting with the lowest angle of attack and 
increasing towards the highest angle of attack. ү 3 
Symmetrical airfoil option is used, the lowest angle of 
attack is zero. ТЕ a non-symmetrical airfoil is used, the 
lowest value for the first two Mach numbers should be -180 
degrees and for all the remaining Mach numbers the value of 
LOMACH entered or -30 degrees. 

All the remaining cards for this particular Mach number 
will contain the data pairs. These cards contain ten £ields 
een, the first field Consisiting of nine cclumns and the 
remaining nine fields consist of seven columns beginning in 
cclumn ten and following the same format as the first card. 
The format 1s (F9.0, 9F7.0). This card is repeated as often 
as needed. Columns 73-80 are not read and may be used for 
ccmments. 

This procedure is repeated until all the CL's are 
entered. Once completed the program is ready to enter the 
values for CD. The whole procedure is repeated аааіп, 
starting with the value for NZ representing the number of 
Mach numbers for which CD's will be entered. ٣٢6 than 
one airfoil data deck is to be used, the abcve ргоседиге 
will start over again by reading in the value cr NZ for the 
number of Mach numbers to be entered for values of CL for 
the second airfoil. The number of data pairs for each Mach 
number must be between 2 and 48. The number of "twice the 
data pairs plus two" must be between 6 and 98. 

Prior to entering the airfoil data, the user should 


review tha following input variables. 
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0 ھ٤[‎ (107) - Symmetrical and nonsymmetrical 
am o даса Input ccntrol. 

2. HIALFA (156) - The highest angle of attack for 
which values of CL and CD will be found at all 
Mach numbers. 

3. LOAFLA (157) - The lowest angle of attack for 
which Values Of CL and CD will be found at all 
Mach numbers. 

4. SPAR (103) - Number of segments using spar 
airfoil data. 

5. TIPSWP (158) = Amount of tip sweep in degrees. 

6. TPSWST (159) - Blade segment number at which 
the tip sweep begins. 

ГЕШИВЗРО 30020) = input control variable for spar 
data. 

85 RB(L) (161-175) - "Controls the blade seqment 


airfoil data assignment. 


DESEAR DATA DECK REQUIREMENTS 


The format for spar data are similar to that of the 
airfoil data with the exception that only one set of spar 
data can be entered into the program. Before the program 
will read spar data, input variable number 120, BSPL, must 
have a non-zero value. In addition, variable 103, SPAR, 
must indicate the number of blade segments which are using 
the spar data. The program automatically assumes that the 
first segment is a spar segment. This is further explained 
in the Case Input section. The spar data are the last tc be 
entered into the program. This is an optional input and is 
moc required. If spar data are not inputed, the progran 
Will assume that the one automatic spar section creates no 
lift and has the drag characteristics cf the first airfoil 


section entered. 
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The format for inputting spar data are as follows. The 
Darst card contains the variable WRSPAR in  Е10.0 format. 
A non-zero value of WRSPAR causes an echo printout of spar 
data to occur. The remaining spar data are handled the 
exact same way as the airfoil section data, starting with 
the variable NZ. There is no input similar to that of the 
airfoil section stating how many different spar data decks 
are being entered since only one is allowed. As before, a 
Minimum of two Mach numbers are required to be entered. The 
blade segment printout indicates spar segments by the use of 
О" for that segment. 

Input variables associated with SPAR data are as 
follows. 

1. SPAR (103) - The number of Elade segments using 

Spar data. 
2. BSPL (120) - Input variable which controls the 


input and use of spar data. 


E. SAMPLE BLADE INPUT 


The next several pages illustrates the sample fcrmat for 
a blade which has the following characteristics: (1) an 
echo printout is not required, (2) there are two airfoil 
decks to be read in and (3) the first airfoil deck has nine 


values of Mach numbers for which CL's are to be entered. 
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fee INPUT CASE LISTING 


Ош 
. [Ж 


ч оу (በ £ UU N ጨመ еН 


222 
235-37 
= 
29067 
68 

69 

70 

71 

(17: 

73 


74 
75 
76 
77 
78 


79 
90282 
8 3 


84 
85 


GREP CASE 99٣ 


DESCRIPTION 
Tip Speed 
Radius 

Speed cf Sound 
Air Density 
No. of Blades 
Forward Speed 


5135 84895. (9273) 
Delta X 

Local Twist 
Local Mass Density 
спора 


В 


Local 
Delta 


рва вов Limit 


Initial Beta 
Initial Beta * 
Initial Beta ** 


Lift and Drag. 
Iteration Limit 


Required Lift 
Required Drag 
Lift Tolerance 
Drag Tolerance 


Birsesloment about 
irap Hinge 8 


Second Moment about 
Shart Orientation 


Pitch-Flap Coupling 
Angle (Delta a 


Dead +0 ті 1 
lat Cyclic -Pitch 


48. 


PROGRAM 
VARIABLE 


OMEGAR 
R 

SPS D 
RHO 
XNB 
VEL 

ER 


рх (15) 

ти (15) 
XMASS (15) 
E 

DPS 1 
ЕТЕТ, 

BIN 

BPIN 
BPPIN 

8 0 


RL 

RD 
XETOL 
XDTOL 
Р4ОМ 


SHOM 
AG BGL; CG 
TDL 


DELD 
А 15 


DIMENSION REQUIRED 
FPS YES 
FT YES 
FPS YES 
SLG/CUFT YES 

- YES 
KTS 89, 90 
- YES 

- YES 
DEG 92 
Sever 78,79 
РТ YES 
DEG жж 

= жж 
RAD ж 
ЗАр/5ЕС ጾ* 
RAD/SEC**2** 

= жж 

LB YES 
LB 95 

LB жж 

LB жж 
SLG-FT 38-52 
SLG-SCFT 38-52 
DEG dor 
DEG - 

DEG ха 


— ——— VE 





Z 


ገ. 


тре 





86 
87 
88 
89 
90 
91 
E 
9.3 


94 
25 
96 
uy 


100 


101 
102 
103 
104-105 


106 
107 


108 
109 


110 
111 
2-113 


ШІ. 
115 
116 


11.7 
118 
19 
120 


Longe yclic Pitch 
collective Pitch 
Inflow Ratio 
Advance Ratio 

VEL Control 
Iteration Output 
Linear Twist 


No. of Blade 
segments 


Climb Rate 
Flat Plate Area 
0 option 


Ела pon MS Re-Use 
ica 


MU Iteration Е 


Tolerance 

Beta Tolerance 
Beta* Tolerance 
Spar Segments 


Second Harmonic 
Control Inputs 


Solrdity 


Symmetric 
АГ Р ОГТ! 


Spring Constant 


Damping Constant 


Patter Option 
Shaft Angle 


Desired AI and 
БИ а прашта 


Al, B1 Tolerance 
Tangent Delta 3 


Phase Angle for 
Delta 


Induced Velocities 
Induced Velocities 
Not Used 

ከ11)! Spar Data 


815 
ТЭЭ 
LAMBDA 
МОТ, 
UIN 
SKIPIN 
TWIST 
XNSEG 


PC P D M 
FPAREA 
TOP 
PCNV 


ADIT 


BTOL 

SP TOT 
SPAR 
ASE ZO 


RSL 
SYN 


SEH 
ЕВИР 


ALOPT 
ALL 


RAIS 
RBIS 


ТОБАВ 
TD3B 
0٤ 


LAML 
UVL 


ES. 


DEG 
DEG 


DEG 


БЕИ 
PITS? 


۶6 


FT-L3/ 
RAD/SEC 


DEG 
DEG 


DEG 


DEG 


жж 
жж 
жж 
6, 90 
6, 89 
жж 
23-37 
жж 


75 


жж 








121 


122 


123 


| 224-125 
1226 
27-137 
6 


139-141 


142 
143 
144 
145 
146 


147 
Щоб 
157 
158 
159 
160 
ИКО 95 


Azımurhal Printout 
Indicator 


Minimum Lift 
Curve Slope 


Iteration Gain 
Factor 


Not Used 
Pre-Coning Angle 
Not Used 


Hub Moment Inplane 
Aero Forces 


АСга гк лч, Roll 
and Pitch Angular 
Velocities 


CG Station 


Rotor Center Staticn 


CG Waterline 
Rotor Waterline 


Aircraft Lateral 
Velocity 


Spar Symmetry 
Airfoil Tables 
Airfoil Tables 
Tip Sweep 
Sweep Station 
Rotor Moments 
Rotor Blade 
Airfoil Data 
ک> ک۸‎ 2101108015 


** Program has 


PESI 
ATEST 


IGC 


PCR 


INPL 
PSIS, PHIS 
TIRS 
FSCG 
TSHR 
WLCG 


WLAR 
V SL Y 


ENIMS ER 
HIALFA 
LOALFA 
11. 251 
ТРУБ Т 
TREIN 
RD 


automatic default 


DEG 


RAD 


RAD/SEC 


INCHES 
INCHES 
INCHES 
INCHES 
KIS 


DEG 
DEG 
DEG 


value. 


жж 


ж ж 


жж 


* x< 





ሰጠው INPUT DEFAULT VALUES 


3 
NO 


68 
69 
73 
76 
77 
81 
85 
86 
87 
88 
91 
93 
97 
101 
102 
114 
121 
П 
123 
156 
157 
159 
167 5 


PROGRAM VARIABLE 
DRST 
FIRE 
Er LEN 
XLTOL 
XDTOL 
BGL 

А 15 
8415 
T75 
LAM3DA 
5 ВЕН 
XNSEG 
PCNV 
BTOL 
BPTOL 
TOLAB 
FEST 

A TEST 
GE 
HIALFA 
LOALFA 
TPESNST 
RB 


Sa. 1f ABIT, item 100, 
xx -50 is TOP, item 96, 


uU 


DEFAULT VALUE 
15.0 
95-0 
15.0 
100. 
50.0 
9080 
=1 
Us 
5 
202 
1.0 
1210 
1.0 

. 001 

. 001 
0225 
DPST 
5.0 Же 
1.0 
30. 
11209 
165 
1840 


LS non-zero 


15 non-zero 


coer 


~ з 


c 
= 
u 
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ШІ CASE INPUT DATA 


Many of the case inputs are self-explanatcry by their 
name listing alone, however, many are not. This section 
will explain the input variables and case options available 
to the user. 

The program enters all the variables into a 200 element 
array called V(I). Prior to entering case data, the program 
will automatically do two things. аи initialize the 
array V(I) to a value of zero. It will assign the default 
values listed in the previous section to those particular 
variables. 

The V(I) array is associated with the variable names by 
equivalent statements. The user needs only to enter values 
for the variables that are different from the default or 
initialized values. In a multiple case computer run, the 
user need only enter variables that are different from the 
preceeding case. If no new value is enter for a variable, 


the value from the previous case is carried over. 
Velocity -_VEL 


The computer program will accept forward velocity in 
one of two ways. The user will input either VEL, item 6, in 
knots or advance ratio MUY, item 89. The value of UIN, iten 
90, determines which variable will be used. If UIN is zero, 
VEL will be used. If UIN is non-zero, MUL will re used. If 
VEL is used, the program calculates the advance ratio by the 


following expression. 


ЯГ" 
حر‎ = 5 (ху сао | A 


If MUL is used, the program will calculate the flight 


velocity by the following expression. 


сл 
(ኣን 
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The offset ratio of the flap hinge, E/R, is the 
ከ r arco from the center of the rotor shaft tc the vertical 


Sipping hinge, normalized by the rotor radius, item 2. 


free becom Seo 2292 - рх 


Delta X is the non-dimensionalized width of each 
individual blade segment starting with segment number one. 
There may be up to 15 segments entered. The number of 
widths. entered here must equal the value of Item 93, XNSEG. 
XNSEG is the number of segments  intc which the blade is 
divided. This can range from two to fifteen. ІЛЕ 15 
recommended that a value of ten or more be used for XNSEG. 
If XNSEG were equal to 12, values cf Delta X would be 
entered for items 8 to 19 and no values would be entered for 
items 20 to 22. The sum of ER plus the summation of the 
Delta X's must equal one. ER is the non-dimensional width 
between the rotor shaft and the flapping hinge cffset. Item 
number 8, which is the first segment width, represents the 
Width between the flapping hinge offset and the point where 
the actual rotor blade airfoil begins. This area is known 
as the spar or cut out segment if no 5раг data were entered. 
If item 103, SPAR, is zero, the program will assume that 
this first section creates no lift and has the drag 
Characteristics of the first inputed airfoil data section. 
Since this area experiences relatively low dynamic pressure, 
the calculations in this segment dc nct have an appreciable 
effect on the outcome of the program. The area between the 
Shaft and the hinge offset, BR, produces neither lift nor 


drag in the program's calculations. The last segment is 


S. 


considered a tip 
to be zero, 
previous runs, the width of this 
to three 
ЖЕР 


percent of the rotor 


The progran 


НЕ into the calculations. If 


item 92, TWIST, can be used. 


the twist can be entered in items 23-37, 


entered for item 92, the 


and ignore 


loss factor segment. 


but drag is calculated in the normal manner. 


Да 


has two options fcr 


If the twist 


program will assume 


all values entered for 


lift is assumed 
In 


The 


section has been 
or Delta X 


set egual 


adius, equal to 


entering geometric 


the blade has linear twist, 
15 non-linear, 
TW. If a number is 

linear twist 
The local at 


TW. twist 


the center of each segment can be entered starting with ¡ten 


MER” If the 


be entered and no values for items 33-37 would be 


A word of caution 1s 
928 


percent chord 


twist option, item 


at the 75 edt te 


assuming that the twist starts at the rotor shaft 


Mencarly out to the rotor tip. 


airfoil started at the 25 percent 


value of -9 degrees, a value 


degrees would 


blade contains ten segments, 


necessary regarding 


The twist 1S considered 


of linear twist equal 


have to be entered 


items 23-32 would 
entered. 

the linear 
to be zero 
The twist is calculated 
and varies 
If the actual rotor blade 
radius pornciyith a tyist 
2 
the 


inp cider for correet 


msc Gistribution to be calculated by the progran. 


5. Items 38-52 - XMASS 


CUHBERED — — ee 120010) 02007 “SEE amp < о — 


The program 


local mass density Or 
The individual mass density for 
38-52, 


the 


in items 
SMOM, 
and 79. 
about 


OL 

about 
тт 

the 


Flapping Hinge, 
data are entered 
Flapping Hinge, 


۰۱ 6٣ ٣ 


94, 


provides two 


the moment 
tha First and 
for item 79, 


the 
provided by items 78 and 


methods fcr entering the 


Of inertia Information. 
each section can be entered 
FMOM and 


items 78 


Second Moment, 


Gan be әпсетей 1 
the Second Moment 
мъ 1 


variable 79 


prcgram use the 


79. ТЕ is 


AA AA 


ና Do cif ае 0 
100 a ar g 

есу! woke 

IS am за 143 

| 





-515"25 0815. 87Н5Л1 
I zp eux Cam 


ШІ то zero, the program will calculate the First and 


Second Moments for items 38-52 and ignore any value entered 
for items 78 and 79. 


6. Items 53-57 - Chord 


The program has two methods for entering local chord 
data. The local chord at each segment can be entered 
Starting with segment number one. If the chord is a 
constant chord, the amount of data to be entered can be 
reduced to only items 53 and 54. If items 53 and 54 аге 
equal, the program assumes that the chord 15 constant 
puroughout the radius and will ignore any “information 
۲١۰٠ ٠ء‎ in items 55-67. Therefore, for constant chord, only 
enter values for items 53 and 54. If the rotor solidity is 
not inputed in item 106, the program will compute solidity 


as follows. 





7. Item 68 - DPSI 


DPSI, or Delta PSI, 1s the incremental azimuthal 
value by which the program advances in its blade flapping 
and force summation routine. This nunter must divide evenly 
into 360 degrees. A minimum value of five degrees is 
permitted. It has been found that for most cases decresing 
the value of DPSI below 15 degrees does not improve the 
accuracy but does increase the computational time. As an 
example, the rotor horsepower required for one particular 
run was 1096 RHP for DPSI equal to 15 degrees, 1095 RHP for 
DPSI equal to 10 degrees and 1097 RHE for DPSI equal to 5 


degrees. DPSI has a default value of 15 degrees. 


0,7 


r e a, 
ag 
45523 








5109 


“318 


FTRL is the maximum limit on the number of times 
that the program will enter the flapping iteration routine 
in.search of a steady state flapping sclution. The progran 
has an automatic default value of 15 iterations. The 
program usually arrives at a steady state flapping solution 
within three to four iterations. The only method by which 
the user can actually determine the number of  flapring 
iterations required would be to use one of the debug output 
options. However, these options will create a huge ancunt 


of output and the user is cautioned abcut their use. 


ከጠ tems 70272, BIN-BPPIN 


— — —— eg ee 


The program has three opticns on how to assign the 
values for the inital flapping angle, velocity, and 
acceleration at the PSI equal zero position. The user may 
(1) input the values, (2) have the program itself calculate 
initial values, or (3) accept the default values of zero. 
Option three is the most used option. It has been 
discovered that there is little or no difference in solution 
time between option two and three. Variable nunber 97, 
PCNV, controls which option is to he used. ТЕ PENV 15 
non-zero, the program will use either option one or three. 
The program initially assigns values of zero to BIN, BEIN, 
and BPPIN before the initial case data are entered.  PCNV is 
assigned the default value of опе. Therefore, with no 
values entered by the user for itess 70-72 and 97, the 
program will start with an initial value for flapping angle, 
velocity and acceleration at the PSI equal zero position of 
Zero. 

If the user sets PCNV equal to zero, the following 
formulas are used for determining the inital values. These 


formulas are derived from page 194 of Ref. 4. 











"муру 





© а > _ AC3e- 221 
в = аа == 





— 
2. 
e E А 
А Se М 
Ser 
.. Ja C 8l» e - 2A] 
 — | 


The term Y is referred to as the Lock number where 


% = C as Qe ጺን 
1% 


In a run where more than one case 1S executed ata 
time, a non-zero value for PCNV allcws the previous case 
values for flapping angle, velocity and acceleration tc be 
used as the initial estimate for these variables in the next 
case. If PCNV is equal to zero, the inital values will be 
calculated by the above fcrmulas. Most cases are run by 


accepting the default value of one for PCNV. 


utem 73 - ХІТТІМ 


= aS ee 


After the program calculates a steady state flapping 
solution for its estimated values of inflow ratic, pitch 
angle and shaft tilt angle, the program determines the 
forces and moments generated by this solution. ЇЇ the 
forces and, optionally, moments do not meet the required 
amounts entered by the user, the program will calculate new 
values to re-enter the flapping routine.  XITLIM determines 
the maximum number of times the program will compare the 
calculated values to the desired values. The program has a 
default value of 15 for XITLIM. A majority ої the 
soluticns require approximately five iterations to converge. 
Once the XITLIM limit is exceeded, the program will stop and 
printout an "Exceeded Limit" statement. The program 
automatically prints as output the number СЕ Mäjor 


Tterations it uses in calculating the solution. 


C 








ПЕ Ета ООШ controls the Short Iteration Option. 
ТЕ XITLIM is a negative number, the program uses this 
option. A complete describtion can be found in the section 


entitled Main Iteration Opticns. 


mm am i ee — COLO ww ee ee ےس ہا‎ J -- 


The required lift, RL, should equal the actual 
weicht of the helicopter that the rotor system is 


Supporting. 


The required drag can be entered in one of two ways. 
It can ре entered directly in item 75 or it can Бе 
calculated by the program by entering the aircraft's flat 
plate area in item 95, FPAREA. The drag is equal to the 
flight path force that the rotor system must overcome to 
Sustain level flight at a certain velocity. Items 75 and 92 
can be entered either aS positive or negative numbers and 
the program w11ll provide the correct fcrward flight solution 
by assigning the proper sign value internally. If a value 
is entered for FPAREA, the program will ignore any value 
assigned to RD. The drag for FPARZA is calculated as 
Follows. 


Б; = AREA 
т) 28 zi ем (FOAR > 


If item 92 is not entered, the user must supply the 
value for the required drag by using the above fcraula. 
There are certain options for which the program igncres or 
does not iterate on drag. An example of this would be when 
the program is used to estimate wind-tunnel test results. 
If these options are desired, see variables 96, TOP, and 
ДО АТОРТ. 
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эг 200101 are the lift and drag tolerances. 
The program will iterate until both the lift and drag are 
Within the tolerances given by variables 76 and 77. The 
program has automatic default values cf plus and minus 100 


pounds for lift and 50 pounds for drag. 


14. Items 78-79 - Moments 


Information regarding the use of FMOM and  SMOM can 


be located in paragraph five on XMASS. 


15. Items 80-82 - Shaft Axis 


The gravitational or weight vector must be oriented 
to the shaft axis of the rotor system. Figure 9 shows the 
positions of these angles. The shaft can be oriented in any 
desired direction. The program will automatically assign 
the proper values for normal helicopter flight. The default 
values for AG, BGL, and CG are 0, 90, and O degrees, 
respectively. This orients the shaft in the vertical 
direction for normal flight. If any other orientation is 
desired, the user must enter the appropiate values for items 
svete 82. 


Delta three inputs are controlled by variables 83, 
115, and 116. These are TD L, TD3B, and PEED, 
respectively. If the flapping hinge is connected in such a 
manner as to cause the blade to change pitch due to 
flapping, this is referred to as a Delta Three hinge. Item 
Ae thE pitch flap coupling angle; 115 is the tangent of 
the Delta Three Bar; and 116 is the phase angle for the 
Delta Three Bar. The program reduces the pitch angle ata 


particular azimuth and segment by the quantity TD3 where 


d. 











Figure 9 Resolution of Gravitational Force 


Control Ants = 
Same Pos - VS 


12 21203 
о ем С 
M Contraco Ахл< = ад 






Tie Parma PLANE 


Control МУ" 5 
SHAFT Ax\s 
Figure 10 Shaft and Control Axis Diagram 


60. 





TD3 = TD3L + TD3B*sin(PSI + PHD3D) 


Normally, variables 83, 115, and 116 are zero. 


The drag increment is a value of delta CD that is 
added to the value of CD obtained from the airfcil input 
data decks. This is added as a roughness factor that 
naturally occurs on blades that are used on production 
aircraft. A value of .002 is normally used. The value is 


not added to the drag calculations for spar data. 


The lateral ande longitudinal cyclic piron = 
controlled by four variables. The program allows the user 
to enter both first and second harmonics of cyclic input. 
The variables 315 and 815, 85 and 86, contrcl the first 
harmonic inputs. The variakles A2S and B2S, 104 and 105, 
control the second harmonic inputs. The A's correspond to 
the lateral inputs and the B's correspond to the 
Mongitudinal inputs. The program has automatic default 
ЕЕ ог 115, 815, А25 апа B2S of -1.2, 7.53, 0 and 0 
degrees, respectively. The user may enter different values 
if desired. Unless other options are indicated, the progran 
will keep these cyclic values constant throughcut the run 
Би vary inflow ratio and pitch angle to obtain a final 
solution. The program will automatically calculate the 
postion of the shaft axis by the momentum theory and will 
use the value of BIS to determine the position of the 
Control axis. Figure 10 demostrates this fact. 

There is a different option available which allows 
the program to seek specific values of longitudinal and 
lateral flapping. This requires the use of variables 100 
and 112 - 114. ТЕ this option is taken, the values or AIS 
and B1S are set equal to zero initially and will be changed 


BY the computer in its itérátion of the required flapfihg 
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angles. The program also has an cption to remove hub 
rolling and pitching moments. Once a solution is obtained 
that meets the lift and drag tolerances, the values of 415 
and B1S are varied to reduce the moments. This option is 
controlled by variable 160. Normally, runs ars made with 
the user not inputing values for itens 85, 86, 1089522105, 
1007 12, 113, 114 and 160. 


Variable 87 is the initial value of the collective 
pitch at the 75 percent radius station at PSI equal zero 
azimuthal station. The program has a default value of five 
degrees. This value is used only to initiate the progran. 
The program, under the normal run option, will vary this 
value in the process of iterating for a convergent sclution. 
There is an cption where the pitch angle remains fixed as in 


a wind-tunnel test. This is the TOP option, variable 96. 


The variable LAMBDA controls the inital estimate of 
an uniform inflov. Since the equations of the program are 
done in a gyrocopter mode, inflow is negative when air flows 
Gowmethrough the rotor. This is the normal fcrward flight 
mode. The program has a default value of -.02 for LAMBDA. 
The program will iterate on LAMBDA in its normal iteration 
pountine. The XLOPT, variable 110, oftion will hold LAHEDA 


eonstant. 


21. Items 892-90 


HUL 


Information regarding the use ዕደ ሸሀፒ and UIN can be 


mound in paragraph one on Flight Velocity. 


22. Item 91 - SKIPIN 


a Se ШЕ = = — — SS 


Information regarding the use cf SKIPIN can be fcund 


in the section of Case Optional Output Indicators. 
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Information regarding twist can be found in 
paragraph four on local twist. 


24. Item 93 - XNSEG 


Information regarding XNSEG can be found in 
paragraph three on Delta X. 


The program can be made to calculate a complete 
solution for any given rate of climb or descent. Climk or 
descent rate must be entered in units of feet per minute, 
With positive values for climbs and negative values for 
descents. The program assumes a uniforn down or up flow 
across the entire rotor surface equal to the rate of climb 
or descent. This value is added as an incremental 
correction into the calculations of UP and will effect PHI 


and angle of attack. 


Information regarding the use of FPAREA, flat plate 


area, can be found in paragraph twelve on Required Drag. 


The TOP option is one of the wind-tunnel cpticns. 
If TOP is a non-zero number, the program iterates to obtain 
Ше required lift of variable 74, but will ignore the 
required drag of variable 75. Item 87, the collective pitch 
at the 75 percent radius at PSI equal zero, will Ге hela 
constant. Item 88, the inflow ratio, will be varied in the 
major iteration routine. A non-zero value of TOP will 
meet in a Value of -50 for ቷቲ6ቹ 122, ATEST. ATEST is the 


minimum accepable value for the lift curve slope when option 


ба. 
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96 or 110 is executed. If non-zero values for both TOP and 
ШЕРТ аге entered, the program will do the TOP option. 
Shaft angle, item 111, must ke inputed by the user. 


Information concerning the use of PCNV, the Flapping 
Solution Re-Use Indicator, can be found in paragraph nine on 


Engtial Flapping Conditions. 


Information concerning the use of PRENT, the 
program's main output indicator, can be found in the section 


entitled Case Optional Output Indicators. 


Item 99 is the End cf Case signal card. It is the 
last data variable that will be entered for each case. If 
XEND is a negative number, the program will stop after it 
determines a solution for that particular case. However, 
Since an infinite number of cases can be entered for each 
computer run, XEND also tells the program if there are тоге 
cases to go. ТЕ XEND is equal to 2.0, the program will 
assume that the next case will begin by reading in new 
ad toil data. If XEND is any other positive real nunber, 
the program assumes that the next case will use the present 
airfoil data and will enter only case input data and any of 
the options which normally follow the case input дата. Each 
time that the variable XEND is entered, Бе especially 
careful to follow the format for NNUM for this card. NUM 
is the number of inputs per data card. NNUM must be a 
negative number for this XEND card. It is this negative 
sign on NNUM which actually keys the computer to stop 


reading data cards for a particular case. 
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Information regarding the use of ABIT can be found 


in paragraph 40. 


32. Items 101-102 - BTOL 
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In the blade flapping routine, the program searches 
for a steady State flapping solution for the given 
Soma 05ہ‎ Of inflow ratio, pitch angle, and cyclic input. 
The program compares values of flapping angle and velocity 
at the PST equal zero azimuthal positicn on each revolution. 
If at this position, the difference between the n-th and the 
(п - 1)th revolution values for flapping angle and velocity 
is less than BTOL and BETOL, respectively, the progranm 
assumes that it has determined a steady state solution. 
BTOL and BPTOL have default values of 0.000001 radians and 
radians per second, respectively. If other values are 
desired, the user may enter those values for items 101 and 
102. 


The number of blade segments using spar data can be 
miouted in item 103. If ne spar data are available, no 
value for SPAR should be entered. Por this case, the 
program assumes that the first segment is the area between 
the flapping hinge and the point where the actual airfoil 
begins on the rotor blade. This area is also referred tc as 
the "cut out" segment. In this case, the program assumes 
that this cut out area produces zero lift and uses CD 
information from the first airfoil section for drag 
calculations. If a non-zero number is entered for SPAR, 
Spar airfoil data must be available. Case input variable 
120, BSPL, controls the spar input option. ТЕ spar data are 
to be inputed, BSPL should be assigned a non-zero value. Іп 


a multiple case run, where spar data are initially entered, 
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the variable BSPL is set equal to zerc as soon as the spar 
data are entered. If the program did not do this, the user 
would have to enter a zero for BSPL fer the next case if no 
new spar data are to be entered. If ina multiple case run, 
the user decides to enter new spar data, the variable ESPL 


must be assigned a non-zero value for that particular case. 
34. Items 


Information regarding the use of the second harmonic 
(Ој inputs can be found in paragraph 18 on Lateral and 


Poengicidinal Cyclic Inputs. 


Information regarding the use Qi ВЕТ retos 


БЕЛО Еу, сап be found in paragraph six cn local chord. 
36. Item 107 - SYM 


SYM is the non-symmetrical airfoil input control. 
If the user assigns a non-zero value for SYM, the progran 
will assume that all blade airfoil data are non-symmetrical. 
The user must enter values for CL and CD for the complete 
range of angles of attack from -180 to +180 degrees. ፲፻ the 
value of SYM is zero, only tabular values from zero to +180 
degrees need to be entered. The above holds true also for 
variable 147, SYMSPR. SYMSPR applies to the spar data 
exactly in the same manner as SYM applies to the airfoil 
data. 


See tons. |06=109 - SPH FDMP 


Values for the spring constant, SFH, and damping 
constant, FDMP, about the flapping hinge can be entered 1£ 
known. These variables can be entered to simulate a 


hingeless rotor system or a system with flapping sptings. 
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This is one of the wind-tunnel options. If ALOPT is 
a non-zero input, the program will iterate to oktain the 
lift required of variable 74 but will ignore the required 
drag, variable 75. In this option, variable 87, collective 
ENTER, will be varied, but not variable 88, inflow ratio, in 
the program calculations. This is the ovposite of the TOP, 
varialbe 96, option. 

If a lift curve slope less than ATEST, item 122, is 
calculated while using this option, the program will stop 
and produce the following message. "Stall Criterion has 
been violated -- will go to next case, if any." tem 111, 
the shaft angle, must be inputed. If non-zero values are 
entered for both TOP and ALOPT, the program will do the TOP 


Option. 


39. Item 111 - Shaft Anale 
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Item 111, the shaft angle, must be inputed whenever 
the TOP or ALOPT options are used. 


40. Items 112-114 


If variable 100, ABIT, is non-zero, the program will 
iterate the blade flapping solution in an attempt to obtain 
the desired lateral and longitudinal flapping angles 
indicated by variables RA1S and RBIS, respectively. The 
program will iterate to the accuracy indicated by TOLAB, 
ከ ን ከን tis Item 112 is RA1S and item 143 15 881۰ 


41. Items 115-116 - Delta 3 


Items 115 and 116 are the tangent and phase angle of 
a Delta 3 Bar. Information regarding the use of 3 TD3B anc 
PHD3B can be found in paragraph 16 on the Pitch-rlap 
Coupling Angle. 
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The program allows the user tc induce a velocity of 
any form onto the rotor systen. This is done in a harmcna#e 
series of the form (А0 + Al*cos(PSI) + B1*sin (PSI) + 
A2*cos (2*PSI) + B2*sin(2*PSI) + ....) for each segment that 
the blade is divided into. If variable 117, LAML, is a 
non-zero number, the program will enter the harmonics of the 
induced velocities. During the Read routine, LAML, will be 
assigned a value of zero. Therefore, for each case where 
different values for the harmonics are desired, LAML will 
have to be set toa non-zero number. Variable 118, UVL, 
controls the use of the induced velocities. II UNE, 15 
zero, the induced velocities will all te set equal to zero. 

A short example will now ре given on the use of the 
control variables in a multiple case run. Assume that no 
harmonic induced velocites are desired for the first case. 
The user would make по inputs for LAML and UVL. For the 
second case, assume that harmonic induced velocities are 
desired. LAML and UVL would be set equal toa non-zero 
number. The harmonic variables would follow the case input 
data for this particular case. For the third case, it is 
desired that the same induced velocites be used. The user 
would not enter any values for LAML and UVL since (1) LAML 
has been automatically set to zero and hence nc new harmonic 
data will be entered and (2) UVL is still equal to a 
non-zero number. It is desired for case four to use no 
induced velocities. The user now would enter the value of 
zero for UVL and program will zero out all the induced 
velocity variables. 

This paragraph will describe how to use the variable 
шаса velocity option. The inputs Al, ВІ, A2, 82, .... 
are numbers in units of feet per second. Ti the direction 
of the velocity is down, a negative sign is asscciated with 


the A's and B's. A negative sign will decrease the angle of 
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attack for an element originally at a positive angle of 
attack. It will increase the amount of negative angle of 
attack for a blade element originally at a negative angle of 
attack. The values of A's and RB's do not effect the 
uniform inflow ratio, LAMBDA, that the program iterates 
upon. The format for using the variable inflow velocity is 
as follows. 

1. A value for NHARM 15 entered іп ІЗ format. 
NHARM is the maximum degree of the harmonics 
that is to be used. If the maximum degree used 
is A3*cos(3*PSI), then NHARM is three. 

The next three paragraphs are repeated for each 
blade segment. 

2. The value for AO is entered in Е15.6 format. 

The IBM 360-67 at the NPS will accept F15.0 
format. 

Бли Тре Values for At, А2, АЗ, ..., ор to NHARM are 
entered on the next data cards in fornat 
ЭЕШ Ц". 

ИШЕ values for B1, 82, ВЗ, ..., Up to NHARM are 
entered on the next data cards іп format 
5E14.6. 

5. Three cards are required for each blade segment 
whether or not the values are equal tc zero. 

ТЕ the user has 15 blade Segments, a minimum of 


45 data cards are required. 


ЦЗ. Тїет_ 120_-_В5Р1 
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Information concerning the use of 8521 can be found 
in paragraph 23 entitled Number of Spar Airfoil Data 


segments. 








44. Item 121 - PPSI 
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5957-011۶6026 5611775 a delta PSI for printout purposes. 
PPSI nas a default value egual to DPSI, variable 68.  PPSI 
Shall never be a smaller increment than the incremental DPSI 


meedeto Calculate the soluticn. 
45. Item 122 


ATEST is the minimum acceptable value of the lift 
curve slope when option TOP or ALOPT is used. If no value 
is assigned, ATEST has a default value of 5 (1/rad). When 
the TOP option is used it has a automatic value of -50 
(1/rad). The lift curve refers to the increase in rctor 
ШЕ With the tilting back of the tip path plane. A check 


on it for a minimum is for ccnvergence purpose only. 


IGC is the Iteration Gain Control factor. If the 
Major iteration fails to converge, choosing a fractional 
value for IGC can greatly speed convergence. This may be 
especially helpful when parts of the rotor are in stall. 
The amounts by which the convergence algorithm changes the 
independent variables is multiplied ty ICC. Setting iten 
91, SKIPIN, equal to zero may heip the user decide if this 


IGC option might be useful. 


ши песен PCR 


The pre-coning angle in radians may be entered here. 


бе веш 136 = INPL 
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Input 1.0 fer INBL in order to remeve hub gesenkt 
inplane aerodynamic forces fron the calculation Of 


aerodynamic pitch and roll moments about the shaft axis. 
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49. Items 139-141 - реве 


Тагтан ёс PSIS, PHIS and THFS represent aircraft 


angular velocites. The aircraft can be given any angular 
velocity in yaw (PSIS), pitch (PHIS), and roll (THFS) in 
radians per second by the use of these variables. They 


rect the calculations of UP, UT and UR.  Faragraph 50 


contains additional information. 
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Tf the angular velccities, items 139 - 141 are 
entered, they are used in the calculaticns of UP, UT and UR. 
Tf no values for items 142 - 145 are entered, the program 
assumes that the rotor system rotates due to the angular 
velocities about the center of the shaft. If values are 
entered for items 142 - 145, the calculations will assume 
that the entire rotor shaft is rotating about the center of 
gravity. Variables 142 - 145 are FSCG, FSMR, WLCG and WLMR. 
FSCG is the longitudinal CG position while FSMR is the 
longitudinal position of the main rotor. WLCG is the CG 
“waterline station while WLMR is the main rotor's waterline, 
or vertical position. All of these values must be entered 


mMm units of inches. 


SW сат 106 - МЕТТ 


The aircraft's lateral speed in knots can be entered 


here. 


52. Item 147 - SYMSPR 


oe ce ee ee el ee um =:: መ= 


If the spar data are symmetrical, do not enter a 
uadue for SYMSPR. The program wili only use values for CL 
and CD between zero to +180 degrees. If the spar data are 
Bon-symmetrical, enter any non-zero value for SYMSPR. The 


program will require values from -180 to +180 degrees. 
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In order for the program to properly calculate the 
region in and around the reverse flcw region, the values of 
CL and CD are required to be known at high and low angles of 
attacks approaching 180 degrees. However, in order to save 
the user from having to enter a whole range of angles of 
attack for all Mach numbers, the program has an option where 
for angles above HIALFA, 156, and below LOALFA, 157 пе 
Mach number is set equal to zero for table lcokup purposes. 
The user is only required to enter large angles for the 
first two Mach number tables. HIALFA and LOALFA have 


default values of +30 and -30 degress, respectively. 
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The GRP prcgram will properly calculate the airflow 
Sweep angle, UT, UR, and pitch angle of a Swept tip airfoil. 
TIPSWP is the amount of sweep of the tip measured in 
degrees. TPSWST is the blade segment number at which the 
Sweep begins. The program assumes that the remaining 
outboard segments starting with TPSWST are swept the number 
of degrees indicated by TIPSWP. UT and UR are modified for 


these segments as shown is Figure 11. 


55. Item 160 = TRIM 


ТЕ TIM 15 a non-zero value, the program will 
attempt to adjust A1S and B1S in order to reduce the rolling 
and pitching moments to less than plus or minus 100 pounds. 
It will first obtain a solution which will satisfy the 
шесшагеа lift and drag, then it will adjust AIS and BIS to 


reduce the moments. 


ОО ከከ ከ 161-175 = ВЗ 


Variables 161 - 175 assign the proper airfoil data 


to blade segments one through fifteen. The program will 
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Figure 11 Tip Sweep Calculation Diagram 
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Бер UN to five airfoil data decks. Тһе RB array is 
initialized with the value of one for all 15 segments. If 
only one blade airfoil deck is used, there is no need to 
enter any values for RB. Only enter values for the segments 
which will use airfoil data sets two through five. I£ no 
Spar data are entered, SPAR equals zero, segment number one 
is considered the cut out segment and is assigned the value 
of RB(1) equals one for calculation purposes and the value 
ERB (I) equals zero for printout indentification. In this 
Case, segment one produces only drag kut no lift. If SPAR 
is greater than zero, RB(I), I = 1, SEAR, will be assigned 
МЕС valle of zero for printout  indentification. These 
segments will use the spar data entered and will calculate 
ШЕН lift and drag. 


E CASE INPUT FORMAT 


Below is located the format that is used to input case 
шаға го the program. All input data cards are cf the format 
IM IJ, OF 12.0) . The input cards contain seven fields 
which are called NNUM, NLOC, С(1), C(2), C(3), C(4) and 
C(5). An example of a data input card is as follow. 

> 1 700. 2628 1148.6 -002246 4.0 

1. NNUM is the number of inputs on this card, 

where C(N) are the inputs.  NNUM must appear in 
either cclumn one or two.  NNUM has a minimum 
of one and a maximum of five. In the above 
example NNUM is five. 

2. NLOC is the item or variacle number of input 
С (1). This refers to the item numbers that are 
found in the section on Case Input Listicgs. 
NDOSEEIS in ІШ format and must be fright 
ОЧЕ зро іп columns four -through six. In the 
example  NLOC refers to item number one, the 


Rotor Tip Speed. 
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ЭЭГ through (. (2) are the values corresponding 
tothe input itens NEOC through NLOC  * NNUM. 
Each value must contain a decimal point and be 
consi 18 for C(1), 19 - 30 for C(2), 
31 - 42 for C(3), 43 = 54 for C(4) and 55 - 66 
for C(5). In this example, values are entered 
for variables one through five, the Tip Speed, 
Radius, Speed of Sound, Air Density and Number 
of Blades. 

4. Omision of NNUM will cause prcgram termination 
with an error explanation statement. NNUM 
greater than five will cause unknown problems. 
Omission of NLOC will cause the present card's 
values of C(N) to be entered into the items 
indicated by the previous card's NLOC. Failure 
о ብብ. 1ህ።ሪፕ ከነ ENLOC, Or NLOG’ greater than 
200, will cause unknown prcblens. Failure to 
properly locate correctly any input value 
Within its own field on the card will cause 
errors in both that input and the number whose 


ле ја encroaches on. 


J. CASE OPTIONAL OUTPUT INDICATORS 


The program has two variables that control the output 
Berntout, varıable 91, SKIPIN, and 98, PRINT. The progran 
will automatically produce a printout cf the case input data 
for each case and a one-page summary of the inital 
conditions and iteration limitations the user placed осоп 
the program. It will also produce a one-page summary of the 
resulting forces, moments and calculated rotor horsepower 
for the final converged solution. The user can also receive 
echo printout of thewairfoil ard spar data decks. If 
Hus echo printout is desired, the user will find the 


Ccrrect printout indicators described in sections C and D 
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entitled Blade and Spar Data Deck Requirements. 

The main optional printout variable is variable 98, 
PRINT. PRINT can be a number from 1 to 1,111,111 depending 
upon the option desired. If no value is enter for PRINT, 
the user will receive the printout described above. Below, 
is listed the PRINT Options. If, fcr example, PRINT is 
assigned a value of 111, the user will receive printout 


options 1, 10 and 100. 
OPTIONAL OUTPUT INDICATORS 


Option 1 Angle of attack, Mach number, 
section Lift and draq 
coefficients, inflow angle, lift, 
and sweep angle at each azimuthal 
position for each radial blade 
segment. Only for converged 
flapping solution.: 

Option 10 Converged flapping angle, rate, 
and acceleration at each 
azimuthal positicn. 

Option 100 Converged integrated forces on 
blade at each azimuthal station. 

Option 1000 Harmonic analysis of blade forces 
for converged case. 

Option 10000 Harmonic analysis of air loads 
for converged case. 

DEBUGGING OR TRANSIZNT OPTIONS | 

Option 100000 Transient flapping angle, rate, 
and acceleration at each 
azimuthal station. 

Option 1000000 Option 100000 plus blade 
velocities, angles, Mach number, 
section coefficients, and lift 
for each blade segment at each 


azinw@enal station. 
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The user is cautioned that the debugging options 


can give a huge amount of output data. 


The second printout option is the variable SKIPIN, 
number 91. This variable controls the summary fcrce, mouent 
and horsepower output discussed in the first paragraph. ፲፻ 
SKIPIN is greater than zero, this summary will be printed 
ВИТУ for the final converged solution. If  SKIPIN is equal 
to zero or a negative number, this summary will be printed 
for each loop through the major iteration (force summation) 
routine. SKIPIN has a default value of one. If the progran 
is not converging to a solution, the user can see 
immediately, with very little extra printout, exactly what 
intermediate solutions the program is producing by setting 
SKIPIN equal to zero. This may help the user in deciding 
whether or not to use the [Iteration Gain Factor, IGC, 


variable 123. 


Нягт 360 EXECUTION CONTROL CARDS 


This section illustrates the control cards required to 
execute the GRP program using the IBM 360 at the Naval 
Postgraduate School. The program may be run under OS or 
CP/CMS. There are two ways of running the program under OS. 
The first is to run the entire program and data through the 
computer at the same time. The seccnd way is to store the 
Main program on a disk as a library program and enter only 
the data through the card reader for each desired case. The 
second method has the two advantages of (1) not requiring 
the user to enter the entire 1100 card main prcgram throuah 
the card reader for each run and (2) the amount of CPU time 
required can be reduced since the main program does not have 
to be recompiled for each run. The program requires 
approximately one minute and forty seconds of CPU time to 


compile. The normal run time for =ach case is approximately 


r r 


nah ог 20: 7 
i =, “пигза NOPY г зуЇр'ї 










Е | ев edt | 
a 21062 
“en al’ 
та са 
i 9 2; 607 
|” ۱ x 
| 9v 12 
є 4% 
> | ú um ЕЧ 
E | « 2 «wt 
е 
= l ru. а 
ӨРТЕ 1 4,72 
» 34 . 34 u ? 


° u A а: 
26-22) 453 "оз 
IA ESA 


IA E e900 LSA 
222000410 мі 77: OM: Ire 0S 
ine Eo. бэвь 355 дећеез ЋЕ 


Г 


.ጀ* -2, ፤ ፻፲ gd дао lr‏ و 


- 


api ара e ^ лой. Deliqwooez 
ЭРТ =“ “FOS SUF гар 14618 


E | = ü 
ин - 22” “ከ22 агч байаа тае VN 


20 CPU seconds. This can vary with the amcunt of printout 
data requested. If the program 1S compiled on the CP/CMS, 
the user will have to request 344K bytes core size on the 
login message. The standard 256K bytes core size is not 
large enough for compiling. However, once the program has 
been compiled, it can be executed within the 256K normally 
available. 

The following cards are required to execute the entire 
program through the card reader at one time. 

Standard Job Card 

АВЕС TORTELG, REGION. FORT=150K 

// REGION. GO=180K 

M FORT. SISIN DD * 

Main GRP Program 

yA 

КОСО о Ек Бр ж 

Case Input Data 

/* 

The following two programs are used to reserve space and 
load the program onto a disk. 

Standard Job Card 

// EXEC FGM=IEFBR14 

//LOAD DD DSN=S1395.HELO,UNIT=3330 VCL=SER=DISKO1 

81. NEN, KEEP) ТАВВІ-БЕТРГ- 150 ,5ВАСЕ- (СҮТІ, (1,1, 1)) 

/* 

Standarda Job Card 

И  ВБЕРС EORTCL, REGEON. FORT=130K 

//FORT.SYSIN DD * 

Main GRP Progran 

Ди 

Ж ЕГЕП СҮЗІП Рр DD 011Т7-3330,Ү01-5ЕНК-015КО01, 

// DSN=S1395.HELO (GRP) ,DISP=SHR 

The $1395 used above and below must be change to S for 
student or F for faculty with the appropiate user number 


instead of 1395. The following cards must be used to 
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execute the program once stored on a disk. 
Standard Job Card 
//GO EXEC PGM=GRP, REGION=180K 
//STEPLIB DD UNIT=3330,VOL=SER=DISK01,DISP=SHR, 
// DSN=S1395.HELO 


МАТО ОСА DD SY1SOUT=A ,DCB=(RECFU=SFBA LRECL=33,BLKSIZE=3325 


2: 022001 рр ~ 
Case Input Data 
yA 


J. Sample Program Output 


This section describes in detail the output availabie 
from the GRP program. In addition, a sample computer output 
is included with each describtion. The program will print 
up to ten different tables. Seven of these tables are 
optional and are not automatically printed. The ten tables 
are as follows. 

1. Echo Printout of Rotor Blade and Spar Data 

2. Case Input Data Card Listings 

3. Summary of Input Data 

ü. Debugging or Transient Information (Options 

100,000 and 1,000,000) 
. Summary of Forces, Moments and Horsepower 


Converged Flapping Sclution (Oftion 10) 


Converged Integrated Forces (Oftion 100) 


. Harmonic Analysis of Z Force (Option 1000) 


оо 44 г Y 
6 


- Converged Blade Analysis (opticn 1) 

|) ስባ Analysis of Air Loads (Option 10,000) 

The input variable PRINT, item 98, controls the output 
ӘБ І1-еше four and six through ten. Items two, three and 
five are always outputed. Item one is contrelled by the 
first data card on the rotor blade and spar section input 


decks. 
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በ ከ ከመ ርን ከ ከገ ጋ ርን Blade Printout 


The next раде contains a partial sample Echo 
Printout o£ Rotor Blade Input Data. This printout 
illustrates that (1) the printout was requested (WBLADE = 
10.), (2) there are two airfoil decks to be entered, (3) 
there are nine Mach numbers for which CL's are to be read in 
and (4) the remaining portion of the printout is the values 
of Mach numbers, angles of attack and lift coefficients for 


mme First airfoil deck. 
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The next page contains a sample Case Input Data Card 
па. This is one of the automatic printouts. It is ап 


echo printout of the input data cards 
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3. SUmmary of Input Data 
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The next page contains the automatic Summary of 
Input Data printout. The following information can be seen 
on this sample output. 

a. The blade was divided into 15 segments starting 


from the hinge offset and proceeding outward. 


Ds No values were entered for the local mass 
desnsity, input variables 38 - 52. Instead, values for the 
First Moment, M(M) = 85, and the Second Mcment, MOM - 


INERTIA = 1450, about the Flapping Hinge, input items 78 and 
79, were entered. 

с. Тһе X row indicates the calculated centers of 
each segment expressed in a percentage of the distance out 
the rotor blade. 

dr The BLADE DECK row indicates that the first 
blade segment was considered a Spar segment. Segments two 
through five and twelve through fifteen belong to airfoil 
data deck number one, while segments six tnrough eleven 
belong to airfoil data deck number two. 

e. The rest of the information, with one exception, 
1s a summary of the case input data. The exception is the 
tern THRUST FACTOR. This is the value used to 
nondimensionalize all the calculated forces in the progranm. 
The THRUST FACTOR equals one о. Moments are 
nondimensionalized by the THRUST FACTOR times the radius. 

f: The program checks to see if all the blade 
segments, delta X's, plus the distance from the shaft to the 
hinge offset, E/R, add up to one. If, on the printout, 
SUM(DX) + E/R does not equal one, the user has made a 


mistake somewhere with the delta X's or in the E/R number. 
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4. Debugging or Transient Information 
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The next page illustrates the Debugging or Transient 
Cen out. This is PRINT Option 1,000,000. The information 
available includes flapping angles, rates and accelerations 
at each azimuthal position with a radial position display of 
НО U, PHI, —Angle-of-Attack, Mach Number, CL, CD and 
Lift Per Inch produced. 

This option is generally outputed only when the user 
is experiencing unknown difficulties with the GRP progran. 
The program will output all of the above information for 
every revolution and iteration until a Ccnverged solution is 
obtained or the program runs out of allowable computer tine. 
If desired, the variable SKIPIN, item 91, will provide a 
printout of the Forces, Moments and Horsepower Summary after 


each major iteration. 
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5. Forces Summary 





The next printout illustrates the Summary of Forces, 
Moments and Horsepower. This page is automatically outputed 
for the converged solution. A printout of this summary can 
be obtained for each loop through the major iteration 
routine by the use of the input variable SKIPIN, item 91. 
The following information can be observed from this sample 
Brintout. 

a. The cyclic lateral and longitudinal inputs, А15, 
Pao, 125 and B2S are printed at the tcp of the page. This 
example indicates that there were no second harmonic inputs 
for A2S and B2S, which is normal. 

b. THETA .75 is the pitch angle of the blade at the 
75 percent radius station at the PSI equals zero azimuthal 
position. In most options the program will iterate upon the 
values for THETA .75 in its search for a converged solution. 

c. LAMBDA refers to the converged value for the 
uniform inflow ratio. If arate of climb or descent was 
used in the case, that rate divided by the tip speed would 
have to be subtracted from or added to this value of LAMBDA, 
respectively. 

а. MU(X)S and MU(Y)S are the advance ratios in the 
longitudinal and lateral directions with respect to the 
Shaft axis. 

Е т, СО, он, СІ апа Ср аге the calculated cverall 
Soerricients of Thrust, Torque, H Force, Lift and Drag. All 
of these items are nondimensionalized by the value of the 
purust Factor. 

f. Lift and Drag forces are calculated with resrpect 
to the relative wind axis. Thrust and H forces are 
calculated with respect to the control axis. Z forces are 
calcuated with respect to the shaft axis. X and Y forces 
are calculated perpendicular to the shaft axis. 

ar The Equivalent Drag is the total drag fcrce 


created by the fuselage, flat plate area times dynamic 


88, 





pressure, plus the profile drag created by the turning rotor 
systen. | 

h. The Equivalent P. A., Or Equivalent Flat Plate 
area is obtained by dividing the Equivalent Drag by the 
dynamic pressure. 

ща Alpha(S) is the shaft axis orientation whils 
Alpha(c) is the control axis orientation. The program uses 
momentum theory to calculate Alpha(S). Alpha (C) 15 
determined from the following relationship, Alrha(S) = 
Alpha(C) + BIS. 

3 9۰۲۰٣۶۰٠٦٦٦٦٦0٥6. 75ت‎ are the lift vector 
offset as a percentage of the rotor radius frcm the rotor 
shaft to create the program's roiling and pitching moments. 

K. РМ and RN are the calculated aerodynamic pitch 
and roll moments about the shaft axis. The SHEARS Hub Pitch 
and Roll Moments are calculated from summing the 


aerodynamic, inertia and elastic hub restraint moments. 
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6. Fiapping Soluticn 


The next page contains a sample printout of the 
converged flapping solution. The first item tc appear are 
the flapping values at the PSI equal zero azimuthal position 
for each revolution prior to the converged revolution on the 
final iteration. The values for flapping angle and its 
first two derivatives are expressed in radians. 

The second part of this printout is the actual 
flapping values for the converged revoltion. The difference 
between the flapping angles and rates between the PSI equal 
zero and 360 degree position must be less than the 
tolerances entered in variables 101 and 102. The numbers 
here are also in radians. The last item isa Fourier 
coefficient series for the flapping angle and it is 
calculated in degrees. All calculations are done in 


respects to the shaft axis. This is PRINT Option 10. 


E. 
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Woe FORCE Integration 


The next table is a sample Force Integration Output. 
This is PRINT Option 100. The forces are a printout fcr one 
blade only at a particular azimuthal position. CQ, CQL and 
COD are the Coefficients of Torque, Torque due to Lift and 
Torgue due to Drag. CQ is calculated from CQ = CQD = COL. 
CX, CY and CZ are all related to forces in the shaft axis 
reference systen. eo с Ене Coefficient of Pitching 
Moment due to aerodynamic, inertia and hub elastic restraint 
moments about the Shaft Axis and CLHS is the Coefficient of 
Rolling Moment due to these same forces. [In the printout 
the (B) character is the number of tlades, which in this 
example is four. SIG is the solidity of the rotor systen. 
MAX B*COD/SIGMA is ael adem stall indicatcr. TE 15 
calculated by determining the first azimuthal station 
between the PSI equal 180 degree and 360 degree position 
that has a value of CQD greater than the value of СОР at the 


PSI equal 180 degrees azimuthal position. 
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The following is a harmonic analysis of the forces 
in the Z or shaft axis direction. Z force is positive in 
| the downward direction. This 1S PRINT Option 1000. 
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2- Converged Blade Analysis 


The 


summary of 


next four pages contain the sample overall 


events occuring on the blade at each azimuthal 


and radial position. This is PRINT Option 1. Variable 121, 


PPSI, controls the aZimuthal intervals that are printed out. 


The output items are as follows. 


а. 
D: 
са 
д. 
е. 
Ёо 
а. 
h. 


X - Center location of the blade segment 
ALPHA - Angle-of-Attack 

MACH - Local Mach number 

С - Local Coefficient of Lift 

CD = Lecal Coefficient of Drag 

PHI - Local inflow Angle 

L(LB/IN) - Lift produced per inch on segment 


Sweep - Sweep Angle of airflow 


В. 
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The next printout contains the harmonic analysis of 
the lift generated per inch on each rotor blade segment. 
tres 1S PRINT Option 10000. The airloads are computed in 
two harmonic forms, both containing terms out to and 
including the tenth harmonic. The forms are 

Lift per Inch = 10- Alcosw - БісіпР -  A2cos24 - 
Бао а -2000с0-2022-283-1037 1..2-27- А10с0510% - 
B10sin10Y and 

пасе рег “Inch = 20 - Clisin(y+4) - C2sin (24+ &) 

-C3sin(3W+4) ..... - С105іл (104%% ). 

The Ratio column is determined from the coefficients 
in Column C. From this ratio, one can immediately determine 
which airload harmonic is dominant and the relative 


relationship of this to the cther harmonic coefficients. 
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IV. GRP SAMPLE ANALYSIS 

The GRP program was executed using data representing a 
relatively new rotor blade. The results were ccmpared with 
results predicted by the  blade's manufacturer. Tho poter 
blade was an unsymmetrical blade. The rotor radius was 
divided into three sections. Sections one and three were 
made of the same type airfoil. The blade included a sweep 
tip design. Тһе blade and helicopter configuration analyzed 
are typical of a helicopter that could be used by the 0.5. 
Navy ina LAMPS type mission. 

For the analysis, the blade was divided into 15 
segments. The program used the manufacturer's values for 
the First and Second Moment of Inertia about the Flapping 
Hinge, vice local blade mass densities. The GRP assumed 
шириною inflow for all flight velocities . It used a rigid 
blade analysis while the actual blade does have live twist. 
The program was run at five different flight weights, 
J ung from 16,359 to 20,829 pounds. A flat plate area of 
35.8 square feet was used at all speeds. The GRP was run 
for forward flight speeds of 40 to 160 knots at ten knot 
increments. The program was executed at sea level, tropical 
day condition: The manufacturer's predicted Botor 
horsepower was obtained from his Shaft Horsepower versus 
True Airspeed curves and was corrected to rotcr horsepower 
by using the manufacturer's Mechanical Efficiency curves. 

The results of the analysis are shown in the next 
Several tables. Table VI illustrates a comparison of the 
GRP required rotor horsepower divided by the manufacturer's 
required rotor horsepower. The GRP agreed within an average 
of two percent on the entire range frcm 40 tc 160 kncts. 
The GRP agreed within an average of one percent for the 
cruise range between 70 and 140 knots. It can be seen that 


between 40 to 60 knots there is a much larger difference 
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between the two 
inflow in this region is not 
Also, 
the manufacturer's horsepower in the 
the 
and much of the retreating blade is in the stall 


mixed and irreguiar. 


This 


helicopter, 
region. 
fuselage attitude at 


represents 


Шо is 


required horsepowers. 


region of 


expected that 


this high speed, 


felt that the 
uniform as assumed, but highly 
the GRP results 
150 to 160 knot range. 


were less than 


For 


a change 


would increase 


the flat plate area above what was used in the frogran. 


GRP progran's 
with those predicted by the manufacturer. 


VELOCITY 


40 
50 
60 
70 
80 
90 
100 
110 
120 
130 
140 
150 
160 


maximum endurance 


TABLE I 
СВЕЩЕНЕ 


1105 
1006 

960 

95 

90 
1047 
1142 
1273 
1438 
1643 
1593 
2300 
2567 


WEIGHT 


= 16359 LBS 
MANUFACTURER'S RHP 


1152 
1041 

973 

958 

971 
1034 
1109 
1237 
1396 
1589 
1852 
2216 
2745 


velocities agreed 


exactly 


RATIO (GRP/MAN) 


295 

537 

22 
1.00 
1-01 
1.01 
1.03 
12093 
„03 
203 
02 
. 04 
. 94 


-d = u 





VELOCITY 
40 
50 
60 
70 
80 
90 

100 
110 
120 
130 
140 
150 
160 


MATOS T TY 
40 
50 
60 
70 
80 
90 

100 
110 
120 
130 
140 
150 
160 


TABLE ZEIT 
ӘКЕЗИНЕ 


1185 
1075 
1018 
1008 
1028 
1086 
1176 
1303 
1465 
1669 
19 21 
2221 
299868 


TABLE III 
GRP RHP 


1366 
1227 
1148 
1118 
1126 
02 
1255 


1375 


1530 
1731 
1983 
2290 
2665 


WEIGHT = 


1247 
1122 
1038 
1009 
1028 
1078 
193 
1264 
1445 
1595 
1902 
2237 
2791 


WEIGHT = 


11:24 


1449 
1301 
1186 
1142 
1153 
1]... 
1273 
1387 
1557 
1746 
2038 
2440 


ЦРЕВО 


MANUFACTURER'S RHP 


i296 TES 


MANUFACTURER'S RHP 


RATIO (GRP/ MAN) 


95 
526 
m» 
1.00 
1.00 
1.01 
1.02 
1.03 
1.01 
1.02 
1.01 
2. 
293 


RATIO (GRP/MAN) 


- 90 
. 94 
2! 
„ 8 
. 98 
<98 
E 
299 
38 
> 
2217 
. 94 
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we LOCITY 
40 
50 
60 
70 
| 80 
| 90 
| 100 
110 
120 
130 
| 140 
150 
| 160 


VELOCITY 
40 
50 
60 
70 
80 
90 

100 
110 
120 
130 
140 
150 
160 


TABLE ТҮ 
GRP RHP 


1407 
1264 
1178 
1142 
1149 
193 
1274 
192 
1546 
1747 
1999 
2305 
2682 


TABLE V 
GRP RHP 


1532 
1370 
1269 
1222 
1219 
1256 
1330 
1445 
1598 
1798 
2045 
2354 
2747 


WEIGHT = 
MANUFACTURER'S RHP 


WEIGHT 


I 


1496 
1340 
1231 
1173 
1179 
1221 
1299 
1414 
1566 
1773 
2048 
2477 


196582LBS 


= 20829 LBS 
MANUFACTURER'S RHP 


1627 
1457 
1338 
1267 
1258 
1237 
1375 
1477 
1638 
1866 
2157 
2589 


RATIO (GRP/MAN) 
„94 
„94 
. 96 
.97 
„98 
„98 
„98 
„98 
„99 
„99 
. 98 
„93 


RATIO (GRP/MAN) 
. 94 
294 
. 95 
. 96 
.97 
„97 
„97 
„98 
„98 
„96 
„95 
„91 
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TABLE VI RATIO COMPARISON 


WEIGHT RATIO 

МЕБОСТТҮ 16359 17321 19246 19658 20529 AVERAGE 
40 06 295 . 94 . 94 . 94 : 95 
50 227 • 96 . 94 . ЭЦ . 94 259 
60 299 23 2227 -96 a99 227 
70 1,00 1. 00 „98 2 2306 .98 
80 1.01 1200 298 21 17 SS 
90 1.07 он . 98 239 227 . 99 
100 1.03 1.02 299 „98 97 1.00 
110 7.03 1.03 222) 29 .98 1.00 
120 12-03 1.01 . 98 MOS (9 1.00 
130 1203 1.02 „99 299 . 96 1.00 
140 1.02 1.01 237 „98 295 EN, 
150 1.04 „99 . 94 593 „91 6ھ‎ 
160 „94 293 = 5 = . 94 

AVERAGES FOR ENTIRE SPEED RANGE 

AVERAGE 1.00 229.9 zn 297 2b 1:۲ 


AVERAGES FOR CRUISE RANGE 
mo = GO KNOTS 
AVERAGE 17:02 1.01 -99 798 201 279 
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V. CONCLUSIONS 


The logic and theory used in the GRP program was 
investigated and found to be sound. However, there were three 
discrepancies in the Navy's version of the program that did 
require attention. The calculations in the reverse flow section 
on the rotor were incorrect, the calculation of the chord at the 
75 percent radius position was incorrect and the original Trim 
Option for reducing moments would not work. All of the above 
discrepancies were corrected. 

Three desirable features were added to the GRP program. 

First, the ability to analyze a rotor blade composed of more than 

one airfoil type was added. The program will now accept up to 

five different airfoil data input decks for use in analyzing a 

rotor system. Secondly, the program would only calculate performance 
in level flight. The ability to calculate performance in climbs 

and descents has been added. Lastly, the program will now calculate 
the aerodynamics for a swept tip rotor blade design. 

The results of the sample analysis described in Section IV 
indicates that the program does produce highly accurate performance 
predictions. The averaged GRP rotor horsepower was within two 
percent of the manufacturer's data. The results were within one 
percent when compared in the area of normal cruise flight. The 
GRP, in this analysis, assumed uniform inflow, constant flat plate 
area and a rigid rotor blade. It is felt that while complicated, 
computer-time-comsuming procedures can be taken to reduce these 
assumptions, they are not warranted if the GRP is to be used 


strictly as a helicopter performance prediction program. 
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